BCS-BEC model of high-T c superconductivity in layered cuprates with 

unconventional pairing 
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High-Tc superconductivity in layered cuprates is described in a BCS-BEC formalism with linearly- 
dispersive s- and d-wave Cooper pairs moving in quasi-2D finite- width layers about the CuO-2 planes. 
This yields a closed formula for T c determined by the layer width, the Debye frequency, the pairing 
energy, and the in-plane penetration depth. The new formula reasonably reproduces empirical 
values of superconducting T c s for seven different compounds among the LSCO, YBCO, BSCCO 
and TBCCO layered cuprates. 

PACS numbers: 05.30.Fk, 74.20.-z, 74.72.-h 



o 
o 

s : 

a ■ 



i 

-a 
c 

o 
o 



> 
in 

On 
O 

o 

O 



X 



I. INTRODUCTION 

It seems to be well-established that central to 
high-r c superconductivity (HTSC) is the layered two- 
dimensional (2D) structure of copper oxides and that su- 
perconducting pairing occurs mainly on the Cu02 planes. 
However, the precise nature of the pairing is still the sub- 
ject of intense research. Recent experiments based on 
angle-resolved photon emission spectroscopy (ARPES) of 
underdoped cuprates suggest that bound fcrmion Cooper 
pairs (CPs) form already at and below temperatures 
higher than the critical transition temperature T c [J- 
4]. Furthermore, ARPES studies of the electron spectral 
function of optimally doped Bi2212 samples show that 
the magnitude of the isotope effect correlates with the 
superconducting gap thus suggesting a role of lat- 
tice phonons in the superconducting pairing. ARPES 
data also suggest that the energy gap A (a measure of 
the energy needed to break a CP) displays an unconven- 
tional d x 2_ y 2 orbital pairing symmetry, with a functional 
dependence A = Aocos26* where = tan -1 K y /K x is 
the angle between the total or center-of-mass momen- 
tum (CMM) HK — (HK X , HKy) of paired electrons in the 
CuOi plane and the a- (or x-) axis while Ao is the value 
of the superconducting gap at the antinode (0 = 0, 7r/2) 
Q. This behavior is also apparent in studies based on 
electronic Raman scattering [7| and in determinations of 
the in- plane magnetic penetration depth A a b [1, Q . 

Although the majority viewpoint in the high-T c com- 
munity seems to argue for such non-s-wave pairing sym- 
metry there are compelling dissenting views, particularly 
work within the past few years, by Miillcr 10] , Harshman 
pl| . Klemm [l2j], and many others. In particular, Muller 
concludes that "...recent experiments probing the surface 
and bulk of cuprate superconductors [show that their] 
character is d on the surface and substantially s in the 
bulk." This conclusion has been bolstered by muon-spin 
relaxation (/iSR) experiments with YBCO reported and 
interpreted in Ref . [ll[ . Several authors [3, 0] have pro- 
posed that the doping process could modify the electron- 



phonon interaction and the Fermi surface with a con- 
comitant shift from d- to s-type coupling as doping in- 
creases. The strongest evidence for an s-wave order pa- 
rameter in a cuprate is reviewed in Ref .[12] where several 
c-axis twist experiments on BSCCO along with earlier 
c-axis tunneling between BSCCO/Pb junctions are sur- 
veyed. Ref. [1 51] summarizes many of the problems with 
the so-called "phase-sensitive" tests [l(| in YBCO. Addi- 
tionally, predictions made in Ref. [l7| that a vortex in a d- 
wave superconductor would exhibit a measurable density 
of states in a four-fold pattern emanating from the core 
have not been observed [3 in either YBCO or BSCCO. 
On the contrary, the vortex cores appear consistent [l8[ 
with isotropic s- waves. 

Replacing the characteristic phonon-exchange Debye 
temperature Orj = hiOo/k B of around 400^ (with k B 
the Boltzmann constant and ujd the Debye frequency) by 
a characteristic magnon-exchange temperature of around 
lOOOif can lead to a simple model interaction such as the 
BCS one but associated with spin-fluctuation-mediated 
pairing [19]. 

The so-called "Uemura plot" 20] of data from //SR, 
neutron and Raman scattering, and ARPES measure- 
ments exhibit T c vs Fermi temperatures Tp = Ep/k B 
where Ep the Fermi energy and k B the Boltzmann con- 
stant. Empirical T c s of many cuprates straddle a straight 
line parallel to the Uemura-plot line associated with the 
simple BEC formula T B ~ 3.31h 2 n 2 J 3 /m B k B ~ 0.218T F 
corresponding to an ideal gas of bosons of mass m B — 
2m* and number density n B = n s /2 where m* is the 
effective mass and n s the number density of individual 
charge carriers. The parallel line is shifted down from T B 
by a factor 4-5. This has been judged [2l[ as a "funda- 
mental importance of the BEC concept in cuprates." 

Previous theoretical papers on the possible origin of 
HTSC 0, [22II24] proposed that the phenomenon might 
be rooted in a 2D Bose-Einstein condensate (BEC) of 
CPs pre-existing above T r . and coupled through a BCS- 
like phonon mechanism [25j, originally taken as s-wave. 
As apparently first reported by Schrieffer [26| , the Cooper 
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model interaction |27j leads to an approximate linear 
energy-vs-CMM relation ^vphK for excited CPs propa- 
gating in the Fermi sea in 3D. This kind of dispersion re- 
lation is not unique to the Cooper model interaction. For 
example, an attractive interfermion delta potential [28j in 
2D (imagined regularized [2!| to support a single bound 
state of binding energy B2 ^ 0) leads [13] (for a review 
see Ref.JHI) to the general CP dispersion relation £k = 
£ + Cl hK+[l - (2 - 16/tt 2 )E f /B 2 ] h 2 K 2 /<im*+0(K 3 ), 
where c\ = 2vp/n precisely as with the Cooper model in- 
teraction [3l|. Hence, the leading-order linearity is not 
induced by the particular interfermion interaction bind- 
ing the CPs but is a consequence of the Fermi sea with 
vf =/= and in which a CP by definition propagates. 
Only in the vacuum limit v F — > => Ef = \ra*vp — > 
does that general dispersion relation reduce by inspection 
to the expected quadratic form Sk = £0 + h 2 K 2 /4m* 
for a composite object of mass 2m*. For either inter- 
electron interaction model, the linear term is a conse- 
quence of the presence of the Fermi sea. The forma- 
tion of a BEC of CPs in 2D does not violate Hohen- 



berg's theorem [32] as this holds only for quadratically- 
dispersive particles. The predicted 2D BEC temperature 

is T c cx (n 2D ) 1/2 oc (0 D T F ) 1/2 where n 2D is the CP 
number per unit area. This leads to values of T c that 
are too high compared with empirical values. However, 
these schemes provide a correct description of other rele- 
vant physical properties of HTSCs such as a short coher- 
ence length, a type II magnetic behavior, and the tem- 
perature dependence of the electronic heat capacity [3] ■ 
They also lead to excellent fits of the condensate frac- 
tion curves for quasi-2D cuprates just below T c [33|, as 
well as for 3D and even quasi-lD SCs. To go beyond 
the simple s-wave interaction, an /-wave formulation of 
BCS theory was discussed by Schrieffer [2(| himself and 
studied in considerable detail by Anderson and Morel 
[HJ in the weak-coupling limit. This has been success- 
fully employed @, Hf| to describe thermodynamic and 
transport properties of high-T c cuprates. The <i-wave ex- 
tension in strong-coupling Eliashberg theory is reported 
inRefs.fH. 



II. BCS THEORY WITH (-WAVE PAIRS 

Some aspects of the Z-wave BCS theory [1 [H, 

rel- 
evant to our HTSC model follow. Consider a system 
of electron- (or hole-) pairs formed via a two-fermion 
isotropic potential V near the Fermi surface and with 
kinetic energies e k = h 2 (k 2 — k 2 F )/2m* (with Kkp the 
Fermi momentum) taken relative to the Fermi energy. 
The Pauli principle prevents background fermions in elec- 
tron states just below (above) the Fermi surface from 
participating in the interaction. In the absence of ex- 
ternal forces each pair propagates freely within a layer 
of finite width 5 along the z direction and infinite ex- 
tent on the x — y plane so that its total momentum 
KK = (KKu , hK z ) is a constant of motion. By neglecting 
spin-dependent interactions the total spin S is conserved 
too and for a spin singlet S — configuration the Pauli 
principle requires that the orbital wavefunction be of the 
form ^(ri.r^) = exp(«K|| ■ R||) cos(K z z)Q(r), where the 
relative coordinate r = ri — T2, Rii is the horizontal 
projection of the CM coordinate R = (ri + T2)/2, and 
K z = rnr/S (with n integer). The z-dependence of the 
wavefunction ensures that the vertical flux of the elec- 
tron (hole) pair across the layer boundary is null. Since 
the relative-coordinate problem is isotropic then $(r) is 
an eigenfunction of angular momentum with quantum 
numbers I = 0,1,2, The total spin S = singlet 
eigenstates of the system satisfy the Schrodinger equa- 
tion 



(fr -V)*(ri,r 2 )=f*(r ll r 2 ) 



(1) 



where Hq is the free Hamiltonian V the interaction po- 
tential and £ the energy eigenvalue. For a given CMM 
wavevector K, we may expand the wave function as 

^(ri,r2) — exp(iK|| -R||) cos(K z z) exp(ik-r). (2) 

k 

In momentum space |T]) thus becomes 

(£ _ e k+K/2 - e k-K/2) a k = ^2 ^kk'Ok' (3) 



Here we develop a general Z-wave BCS-type theory 
which is then applied in a quasi-2D BEC picture with 
either I = or / = 2 pairing symmetry. In § II the /-wave 
BCS theory within the framework of the present model is 
discussed. In § III we study a quasi-2D BEC of linearly- 
dispersive, massless-like CPs and we evaluate the number 
density. In § IV the areal density n 2D of charge carriers is 
estimated by calculating the magnetic penetration depth 
arising from the CPs. In § V an analytical expression for 
the critical BEC temperature is derived, which is then 
applied in § VI for various superconducting materials in- 
cluding YBCO under different doping levels. Discussion 
and conclusions are given in § VII. 



with Vkk' = < k, — k|V|k', — k' >. Since the interaction 
potential V depends only on r it admits the expansion 

00 / 

y ™ = E E v ^ ik'i)irwir ro (n k o. (4) 

1=0 m=-l 

For small coupling amplitudes V/(|k|, |k'|) the contribu- 
tions of different I spherical harmonics 17 n (f2k) ^ n © can 
with good accuracy be considered relatively independent 
[34[. In that case ([3]) yields an analytical solution by as- 
suming that Vi is separable, i.e., Vi(|k|, |k'|) = V Q (l) f l k f l k * 
so that §5§ becomes 

{£ - ^k+K/2 - ek-K/2) «k - V^f k ]T a k ,f k l (5) 
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where at = cikY" 1 (^k)- Eq.© can be now rewritten as 



yil) f l 



c{l) 

~ e k+K/2 — e k-K/2 



(6) 



where ay f l k * = is a constant. One thus obtains 
a BCS-type integral relation for a CP in the eigenstate 
characterized by (/, m) 



1 = V r 



(0 



\fl\ 2 



c(l) 

k C K ~ e k+K/2 - £k~K/2 



(7) 



Following Schrieffer [2a] we assume that the angular- 
independent I component of the generalized BCS model 
interaction Q is given by 



(8) 



with Vq > for CPs with relative momenta (k, k') lying 
in the neighborhood of the Fermi surface 



k F < |k + K/2|,|k-K/2| < K n 



(9) 



and V^fifC 



otherwise. Here K n 



k 2 D 



with defined in terms of the Debye e nerg y via ftun 
f?k 2 D j2m* . A straightforward analysis [26[ reveals that 
([7]) yields a bound state with energy £® < for arbitrar- 
ily weak coupling so long as the potential is attractive in 
the region (j9|) in fc-space. Then, a bosonic CP can form 
only if the tip of vector k lies within the intersection of 
the two spherical shells defined by © whose center-to- 
center separation is K\ fermions with wave vectors lying 
outside this overlap are unpairable [24j . 

In the quasi-2D limit the fundamental expression ([7]) 
can be evaluated by substituting the summation over k 
by a 2D integration. In addition, for small 6 the only 
term in K z that yields a finite contribution is n — 0. By 
assuming a 2D cylindrical Fermi surface we obtain 
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kdk 



r l4 } ! 
+ 2Hujl 



e k+K/2 + e k-K/2 

- vphK cost 



-(10) 



Mi 

A" 



vphK cost 



where h = k F + (K/2) cos 9, k 2 = k F + k D - (K/2) cos 0, 
and the approximate equality in the second row holds 
up to terms of order (fco/fep) 2 = Qd/Tf- The angular 
brackets denote an average over a 2D cylindrical Fermi 
surface < ... >f~ > (1/2tt) f Q ™ dd. The Fermi average 
then gives the energy spectrum of excited CPs [261 ]: 



aStf' +c 1 hK + 0(K 2 ) 



7(0 



(11) 



where c\ = 2vp/n in 2D and is the binding energy 
of the CP ground state (hK = 0) ^ 

= -2^/[exp(2/AT T/ (i) ) - 1]. (12) 



The dispersion relation (fTTj) is linear in leading order 
rather than quadratic as would be expected in vacuo. 
As a consequence, all excited CPs behave like a gas of 
free massless-like bosons with a common group veloc- 
ity ci = K^dEg- /dK, but a variable energy determined 
by their CMM HK. The dispersion relation ([TT]) implies 
that in order for a CP to remain bound (i.e., < 0) its 
maximum CMM wavenumber must not exceed the value 



\SiP\/ci = K since CPs with K > K have Eg > and 
thus break up [26]. 

Explicit expressions of relevant thermodynamic vari- 
ables and transport coefficients evaluated within the 
weak-coupling limit of the l-w&ve BCS theory have been 
derived in Refs.Jl, H3, In these papers it is shown 
that the average behavior of most of these quantities over 
the cylindrical Fermi surface exhibits small variation due 
to the explicit realization of an I = or I = 2 symme- 
try [f| [35[ . In particular, the temperature-dependent gap 
equation is given by @, [26|, [13] 



rhuj D ; tanh 
1 = Af V / de k 
Jo 



^/e 2 + A(0 2 | 5 (0(k)|2 



(13) 

where /3 = l/k B T, g (0) (k) = 1 for I = 0, and 5 (2) (k) = 
cos(2#) for I = 2. The critical temperature is determined 
from (HH) by the condition AW(T C ) = 0. In the weak- 
coupling limit fcjfl /ksTc ^> 1 it can be calculated an- 
alytically and it follows that T c is independent of the 
I -state [3^: 

2e 7 

k B T c = hu D exp(-l/M V ) £ii l.l3hu D exp(-l/Jsf V ) 

(14) 

with 7 ~ 0.577 • • • Euler's constant. In the zero- 
temperature limit Ag = A^(T = 0) the energy inte- 
gration in (fTB")) leads to the gap relation [§3] 
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2T^haj D exp{-l/JV V {l) ) 



(15) 



where r« = exp[-(| ff « | 2 In | 5 « |) 



For 



o, r(°) 



1, while for / = 2, r< 2 ) = 2exp(-l/2) ~ 1.213, so that 
combining (fT4")l and (fT5"j) we are led to the gap-to- T c ratios 



2A 



(0) 



ksTc 



3.53 



2A 



(2) 



kBT c 



4.28. 



(16) 



For I — one recovers the standard BCS result [25| and 
the somewhat higher value for the I = 2 d-w&ve case. We 
note that the quantity A® /T® has the same functional 
dependence as the zero-temperature gap of the BCS the- 
ory [25j]. Considering that measurements of the energy 
gap for any given cuprate show some scatter about a cen- 
tral value Aq XP [3^] in the following we shall assume that 

A e xp ~ a 2) /r( 2 ) ~ A (0) 



Ar 



On the other hand, the average superfluid density 
p s {T) = A 2 b (0)/A 2 b (T) exhibits a more pronounced an- 



gular momentum dependence. This is given by 

l5 W (k)| 2 



P 



(0 



1-/3 / de k ( 

Vosh 2 (±Py/el + A(') 2 | ff (')(k)P 



2\ I F 



(17) 

In the low-temperature limit (jTTJ) yields for I = an 
exponential T-dependence 



^)(T)~1 



2ttA 



(0) \ V2 



fc R T 



exp(-A^ 0) /fc B T) (18) 



while for I = 2 it gives the linear T-dependence 



(21n2)fc B T 



A 



(2) 



(19) 



Experiments [U, 0] on the temperature variation of the 
magnetic penetration depth X a b(T) are consistent with 
the quasi-linear behavior (|19[) which is a signature of d- 
wave symmetry. However, its asymptotic value X a b{T — > 
0) is independent of I, a result that we apply below. 



III. BOSE-EINSTEIN CONDENSATION 

We assume that charge carriers are an ideal binary 
mixture of non-interacting unpaired fermions plus break- 
able bosonic linearly-dispersive CPs [14 . 122L |23l |26I | . Let 
the fermion number per unit area be rif = w/i + ri/2 
where n/i and n/2 are the number densities of unpairable 
and pairable fermions, respectively. Unpairable fermions 
lie outside the interaction region of (J5J) unlike the pairable 
fermions whose T-dependent density n/2(T) is 

n f2 (T) = 2[ nt D {T) + n 2 ° K < Ko (T) ] + n u f2 (T). (20) 

Here n^ D is the bosonic number density of CPs with 

CMM wavenumber K = 0, n-o<x<x that with < K < 
K , and n J 2 the number density of pairable but unpaired 
fermions. By asserting that in thermal equilibrium these 
kinds of fermions arise precisely from broken CPs [22| we 
identify n u f2 (T) = 2n 2 £ o<K<Kmax (T). On the other hand, 
at T = all pairable fermions should belong to the con- 



densate (Ref.[38j, p. 122) so that n f2 (0) = 2ng B (0) = 
2n 2D where n 2D is the total boson number per unit 
area. The number equation for pairable fermions may 
thus be reexpressed in terms of boson quantities alone, 



namely n 2D = n 2 , (T)+n 2 



,2D 



(T) 



2D 



2D 



2D 

1 0<K<K„ 



(T) 



<K<K 

Thus 



,2-D 

l K <K<K„ 
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d 2 K 



o+ 



(27T) 2 



"iexp^ - 1 



XT) 



(21) 



where (3 = 1/fcsT, fj, the boson chemical potential and 
z = exp/3/x is the fugacity (0 < z < 1). On in- 
troducing (|2"Tj) the energy-shifted boson dispersion re- 



lation £ 



(') 

K 



hc\K for K > the integral can eval- 



uated by changing to the variable x = j5hc\K. Since 
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FIG. 1: Comparison of experimental T c s vs. theoretical pre- 
dictions (|29p as function of zero-temperature inverse penetra- 
tion length A",, 1 for YBCO compounds with different doping 
degrees. Square datapoints are taken from Ref . [37t] . except 
for uppermost square referring to the optimally doped regime 
[3Sl |. Vertical "error bars" represent full widths of ci peaks, 
where ui is the real part of the conductivity a employed in 
Ref. [37l| to determine A~ 
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Ci = 2vp/-K and K max = kpil + kp/kp), the upper inte- 
gration limit x m ax in (|21[) is then be very large, namely 
Xmax — fitwpkp = 2Ep/kBT 3> 1. The last inequality 
is consistent with the maximum empirical value for the 
ratio k B T c /E F < 0.05 reported [H for all SCs includ- 
ing cuprate SCs. Given the rapid convergence of Bose 
integrals the upper integration limit x ma x may safely be 
taken as infinite in (|2Tj) so that the integrals can then be 
evaluated exactly by expanding the integrand in powers 
of zexp(— x) and integrating term by term. The number 
density (f2"Tj) becomes 



2 00 



(k B T) 

2ttH 2 c 2 ^ n 2 

1 71— \ 



E 



(22) 



The critical BEC temperature T c is now determined by 
solving (|22J} for nl D {T c ) = and z(T c ) = 1. One obtains 



Tr 



where C(2) = 7r 2 /6. 



ha 
kn 



2ixn 



2D 



cm 



1/2 



(23) 



IV. CHARGE CARRIER DENSITY 

The areal density of charge carriers was formerly es- 
timated from measurements of the London penetration 
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superconductor 9o (K)° A (mcV)'' X ab (nm) c 8{A) d T^ xp {K) e Tj h (K) {2A a /k B T c ) expU "> (2A /k B T c ) t ' 



(La. 9 25Sr.07 5 )2CuO4 360 6.5 250 4.43 9 36 36.4 4.3 4.14 

YBaaCuaOe.eo 410 15.0 240 2.15 h 59 56.0 5.90 6.09 

YBa 2 Cu 3 06.95 410 15.0 145 2.15 9 93.2 92.6 4.0 3.68 

Bi 2 Sr 2 CaCu208 250 16.0 250 2.24 s 80 72.2 4.64 4.85 

BiaSraCaaCuaOio 260 26.5 252 2.24* 108 109.2 5.7 4.99 

Tl 2 Ba2Ca 2 Cu2 08 260 22.0 221 2.14 4 110 104.1 4.5 4.47 

Tl 2 Ba2Ca2Cu 3 Oio 280 14.0 200 4.30* 125 105.5 3.1 2.96 



TABLE I: Physical parameters of cuprate superconductors and predicted values for T c , and the ratio 2Ao/fcsT c according to 
(|28|) . Debye temperature is 9d = hujo/ks- Parameters taken from from Ref. 38] (see also references cited therein): a) table 
4.1, b) table 6.1, c) table A.l, d) table A. 2, e) table A.l, f) table 6.1, g) estimated from band-structure calculations |39l. [40] . h) 
estimated as <5 = 0.64 d„t, and i) estimated as <5 = 0.68 Ci n t, where Ci„t is the CuC>2 interlayer separation for a given cuprate. 
For YBCO 6d, Ao, and S are assumed the same for different dopings. 



depth Al which is the distance over which an external 
magnetic field decays within the superconductor. For 
superelectrons with a 3D density n s , charge e, and effec- 
tive mass m* , one has the well-known relation 1/A 2 = 
47re 2 n s /m*c 2 . By introducing [20| the average interlayer 
spacing Ci n t between CuOi planes in HTSCs it follows 
that n 2D ~ Ci n tn s . Penetration-depth data spanning a 
wide range of critical temperatures are consistent with 
the phenomenological Uemura relation T c oc 1/A 2 ^ oc 
n 2D /m* H3- 

Within the present model we evaluate the magnetic 
penetration depth due to linearly-dispersive CPs with 
charge 2e, and constrained to move within a thin layer of 
width (5 with a uniform CM speed c\. Thus, we first con- 
sider the expression for the 3D supercurrent of excited 



CPs \14 J, 



,3D 



(2e)ciK with K ee K/K. We now 



introduce the contour integral of the CP wavefunction 
phase within a homogeneous medium and in the pres- 
ence of an external magnetic field B = V x A. The 
following integral along any closed path vanishes 



HK + -A 

c 



dr = 



(24) 



where c is the speed of light in vacuum. By expressing 
K in terms of J s and using Stoke's theorem to evalu- 
ate we get a modified version of London's equation 
J s = — A p A where A p ee 4e 2 cin 3D /hcK . Taking now 
the curl of this modified London equation and introduc- 
ing Ampere's law VxB = (4ir/c)J s , it follows that the 
magnetic induction B satisfies the Hclmholtz equation 
V 2 B = A" 2 B, where 



1 

A2 



(2e) 2 47tci?t 



(25) 



Note that London's result is recovered for quasi-particles 
with density n 3D —> n s /2, momentum hK — > 2m*ci, 
and charge 2e — > e. This expression for A varies be- 
tween its minimum value A = when K = (perfect 
diamagnetism), and its maximum, say Ao, for K = Kq 
(CP breakup). It seems natural to identify Ao with the 
experimentally observed value of the in-plane penetra- 
tion depth at T = 0, namely Aq = X a b(T = 0). Here, 



Kib (0) — + -V^O) i s t ne geometric mean of this 

parameter measured along crystallographic in-plane di- 
rections a and b. As shown in § II, this parameter is in- 
dependent of the explicit value of the angular momentum 
I. By substituting the dispersion relation (|11[) to elimi- 
nate Kq from Ao and imposing the relation n 2D = Sn 3D 
the 2D charge carrier density becomes 



2D 



Mi 



1 



l^c{ X 2 ab 



(26) 



This latter expression can be reformulated by considering 
the relation (jl~5j) and the weak-coupling limit of (fT2"|) . It 

follows that \£^ l) \ = (A^) 2 /2Huj d so that 



2D 



SA 2 



32irc 2 c 2 hidD A 2 b 



where the approximate relation Aq /T® 
fied in S II was used. 



(27) 



Aq as justi- 



V. CRITICAL TEMPERATURE 

The final explicit expression for the critical BEC tem- 
perature T c is now obtained by substituting ([2T]) in 
This leaves 



he 



35 



2nkBe \2frui£> 



1/2 



^ (28) 



which is independent of the CP speed c±. We observe 
that for fixed values of u>d, Ao, and 5, the critical tem- 
perature increases linearly with X~ h . This dependence 
has been observed by Zuev et ai.[37| in experiments in 
underdoped YBCO films with T c s ranging from 6 to 50JT. 
They conclude that, within some noise, their data fall on 
the same curve p s oc A~ b 2 oc T 2 3±0 - 4 , irrespective of an- 
nealing procedure, oxygen content, etc. Thus, by assum- 
ing that except for A a fc the other YBCO parameters are 
approximately independent of the dopin g lev el, we intro- 
duced in OS) the values: 8 D = 410if (H], A = 14.5 
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meV [38|, and S = 2.15 A [39|, |4(| to get the relation 



16.79[(/im)- 1 j£:] 

A a 6 



(29) 



Figure 1 is adapted from Ref.[37| and compares theoret- 
ical predictions (|2"9"|) with experimental data, as well as 
with data pertaining to higher doping regimes. We see 
that (|29[) gives an excellent fit to the experimental data. 
The same functional dependence has been observed in 
single YBCO crystals near the optimally-doped regime 
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More recently, Broun et al. @ found that their 
samples of high-purity single-crystal YBCO followed the 
rule T c cx A^ 1 oc n~J 2 oc (p — p c ) 1 / 2 where the doping 
p is the number of holes per copper atom in the Cit0 2 
planes and p c the minimal doping for superconductivity 
onset. The measured value of the penetration length in 
YBCO crystals is an order of magnitude bigger than in 
thin films 0, |4l|, so that the specific values of T c s de- 
rived from (|28|) are not in such good agreement as in the 
YBCO films. However, one should expect variations of 
parameters such as the energy gap associated to crys- 
tals and film systems. It has been pointed out 37J that 
YBCO films seem to behave more like other cuprates 
such as BiSrCaCuO or LaSrCuO than do YBaCuO crys- 
tals. Furthermore, a different approach 42] based on 
measurements of the lower critical magnetic field H c i(T) 
for highly underdoped YBCO indicates that experimen- 
tal data may be consistently described only by assuming 
T c oc rig 61 , in close agreement with studies mentioned 
above. 

Theoretical values of T c for superconducting 
cuprates with different compositions have been 
also calculated using (|28j) . Here we report on 
these seven layered-cuprate superconducting com- 
pounds: (La. 9 25Sr. 075)2CuO 4 ; YBa 2 Cu 3 6 .6o; 
YBa 2 Cu 3 6 . 9 5; Tl 2 Ba 2 Ca 2 Cu 2 8 ; Tl 2 Ba 2 Ca 2 Cu 3 Oi ; 
Bi 2 Sr 2 Ca 2 Cu30io; and Bi 2 Sr 2 CaCu 2 0s- Characteristic 
parameters for these materials were taken from tables 
compiled in Ref.[38| (see also {43T - I451 ] ) . Concerning 
the layer width 5 no direct experimental data are 
available. We have employed results derived from energy 
band-structure calculations for cuprates. Contour 
plots [3j| HjJ of the charge distribution for La 2 Cu04, 
YBa 2 Cu307, and BiCa 2 SrCu 2 0s suggest that charge 
carriers in each of these systems are concentrated 
within slabs of average width 5 ~ 2.61A, 2A5A, and 
2.28A, respectively, about their Cm0 2 planes. As Ci n t 
denotes the average separation between adjacent Cu0 2 
planes, it follows from crystallographic data [HI that 
the yttrium and bismuth compounds give S ~ 0.64 Ci n t 
and 0.68 Ci n t, respectively. Taking into account that 
BiSr 2 Ca n Cu n +i06+Ti compounds possess the same 
layeringscheme as their TlBa 2 Ca n Cu„+i06+r; counter- 
parts j38|, we assumed that the condition S ~ 0.68ci„ t 



holds also for the thallium compounds. The former esti- 
mations are congruent with Uemura's surmise [20l | that 
SC charge carriers in layered cuprates are concentrated 
within slabs of width 8 = Ci nt . 

Table I shows results obtained using the foregoing as- 
sumptions, together with the physical parameters in- 
volved in the calculation. In most cases we find rather 
satisfactory agreement between predicted and measured 
values of T c . We also find very good agreement 
between theoretical and experimental gap-to- T c ratios 
2Ao/ksT c . Average theoretical and experimental such 
ratios presented in Table I are (2A /fcsT c )*' i ~ 4.45 and 
(2Ao/kBT c ) exp ~ 4.59, respectively. Both are consis- 

(2) 

tent with the ratio 2Aq /fcsT c ~ 4.28 predicted by the 
I = 2 BCS theory in (fTB"|) . We have not attempted es- 
timate uncertainties of our theoretical results since the 
accumulated data of the physical parameters involved in 
the calculation, particularly Ao and A a b, show a wide 
scatter. 



VI. DISCUSSION AND CONCLUSIONS 

We have shown that layered-cuprate HTSC can be de- 
scribed by means of an i-wave BCS theory for a quasi-2D 
BEC of Cooper pairs. The theory involves a linear, as 
opposed to quadratic, dispersion relation in their total 
or CM momenta. The theory yields a simple formula for 

the critical transition temperature T c with a functional 

l li 

relation T c oc 1/A a fc oc n s which applies to a variety of 
cuprate SCs over a wide range of dopings. Although this 
behavior apparently disagrees with the phenomenological 
Uemura relation T r oc 1/A^ h [201 ] . different experimen- 
tal studies 0, [13, EH show consistency with the inverse 
linear dependence of T c . Additional consistency is also 
seen with the reported dependence T c oc n"' 61 arising 
from measurements of the lower critical magnetic field 
[42l ]. When averaged over a cylindrical Fermi surface, 
the physical quantities involved in the theory show small 
dependence on the angular momentum state I. However, 
the gap-to- T c ratio 2A /fcsT' c is closer to that predicted 
by the extended BCS theory for I = 2 than for I — 0. It 
is shown elsewhere [46| that all relevant 2D expressions 
derived here arise in the limit ksTS/hci — > of a more 
general 3D BCS-BEC theory for layered materials. 
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